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Abstract: We study the dimensional asymptotics of the effective actions, or 
functional determinants, for the Dirac operator D and Laplacians A + PR on 
round S"'. For Laplacians the behavior depends on "the coupling strength" /3, 
and one cannot in general expect a finite limit of C'(0), and for the ordinary 
Laplacian, (3 = 0, we prove it to be +00, for odd dimensions. For the Dirac 
operator. Bar and Schopka conjectured a limit of unity for the determinant 
([BS]), i.e. 

hm det(A^c"an) = l- 

n — *oo 

We prove their conjecture rigorously, giving asymptotics, as well as a pat- 
tern of inequalities satisfied by the determinants. The limiting value of unity is 
a virtue of having "enough scalar curvature" and no kernel. Thus for the im- 
portant (conformally covariant) Yamabe operator, /3 = (n — 2)/(4(n — 1)), the 
determinant tends to unity. 

For the ordinary Laplacian it is natural to rescale spheres to unit volume, 
since 

^hm^det(A5--.)-i. 



1. Introduction 



Ever since spectral zeta functions of natural geometric elliptic differential op- 
erators on manifolds started appearing in mathematics ([RSlj, ^RS2j ) and in 
the regularization of path integrals (fHaj), there has been interest in calculating 
values of the associated zeta determinants. Explicit formulae on various spaces 
have been derived ([E], [BS], [DK], [D5^ . [D5^ . [Ki], [^2]), 
and some authors have displayed numerical features of the behavior of such de- 
terminants, culminating in the formulation of Conjecture 1 in [BS', concerning 
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the limit of the determinants of the Dirac operator on standard n-dimensional 
spheres, as the dimension n grows large. Thus it seems an opportune time for 
investigating rigorously the asymptotics, using the previously established formu- 
lae for determinants. Contrary to what one might think, the rather complicated 
looking series, which involve Barnes zeta functions and generalized Bernoulli and 
Stirling numbers, can often be understood via quite elementary methods. 

In Section[2]we find asymptotics of the effective action for the Dirac operator, 
and of the phase factor sometimes included. As a corollary this resolves the 
questions raised in the form of Conjecture 1 in [BSj . The key point in the proof is 
the use of explicit formulas by Thomas Branson, Bar-Schopka and J. S. Dowker, 
and application of recursion to control the indirectly defined special functions 
and polynomial coefficients that appear, along with the functional equation for 
the Riemann zeta function and some elementary estimates of (r ([Br], |BS| . 
[Dowl] , |Dow2| . See also the book by K. Kirsten [Ki]). 

In Section [3] for the Laplace-Beltrami operator on functions, more refined 
methods are needed. The reason seems to be that the ordinary Laplacian is not 
a "natural object" to consider. While an approach using polynomial coefficients 
such as (generalized) Stirling and Bernouilli numbers works well for the more 
natural conformally covariant Yamabe and Dirac operators, where in fact the 
zeta function and derivative at zero tend rapidly to zero, then in the cases 
when it tends to oo, formidable cancellations amongst large terms should occur 
in such sum formulas. Namely for Dirac and Yamabe the convergence holds 
with absolute signs inside the Bernoulli/Stirling expansion sums, while for the 
ordinary Laplacian they do certainly not. As shown below, the zeta derivative 
at zero of the ordinary Laplacian goes to +00 as the logarithm of the dimension, 
while the same quantity for Dirac and Yamabe operators converges exponentially 
to zero. 

The proof in the ordinary Laplacian case instead relies on estimating a broad 
range of Barnes zeta functions through convenient contour integral representa- 
tions. Namely, by exploiting geometry of the situation it is possible to deform 
the basic Hankel contour to make use of a hybrid of Laplace's method and the 
method of stationary phase, the idea of which is exactly that of cancellation of 
increasingly rapid oscillations. We note that the standard formulations of the 
method of stationary phase and Laplace's method for contour integrals (see for 
example [01]) do not seem feasible for the estimates needed here. What is needed 
is a mix of these standard tools, due to phase factors which are neither real nor 
purely imaginary. In this respect the present paper may also turn out useful for 
treating similar problems on contour integrals. 

2. Dirac operator sphere determinants 

As in (BSj . we define the Dirac operator determinant including a phase as follows. 
Definition 1. 

det{D) :=exp(z|(Cz52(0)- 771,(0))) exp(-icD^(0)). (1) 

To state the theorem, we fix the notation that ipo — fCD2(0) — 77_d(0). The 
theorem in particular proves Conjecture 1 of }BSj . 
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Theorem 1. For the standard round spheres we have: 

lim det{D, S") = l. (2) 

In fact we have the following asymptotics 

|det(Am =exp(o((|r)), 

Furthermore we have the inequalities 



detp,^")| > 1, 
\detiD,S")\ > 1, 



(PD > 0, '1 


If n = 


(mod 4), 


ipD ^ 0, 


if n = 


1 (mod 4 ), 


ipD <0, '1 


if n = 


2 (mod 4 ), 


ipD =0, '1 


if n = 


3 (mod 4 )■ 



(3) 



Remark 1. This pattern is also visible in the numerics in [BSj and is interesting 
in comparison to the present author's results on local extremals of determinants, 
where that pattern is found to be (max, max, min, min), again depending on 
the respective dimensions mod 4 (see |M0| ). This is yields another example of 
mod 4 dependencies generically showing up for zeta regularised quantities. 

Proof. We need to show that, writing from now on 5" for S'"ani 

lim C(r>2,S")(0) = 0, 

lim Cm2.s")(0) = 0- 

We firstly review the expressions for the zeta functions of Dirac squared, by 
Tom Branson Br . Depending on the parity of n we have as follows, in perfect 
agreement with 



2fc+i '="1 

C(D^5^")(5) = .^fc - 1)1 X! ^".fcCfi(2s -2a- 1), if n = 2fc, 

\ (^) 

C(D^,S-){s) ^ ea,fc(2^^-"" - l)Cfi(2g - 2a), if n = 2fc-l, 

where the da.k and ea.k are integers defined indirectly through the following 
polynomial expressions in x. 

fe-l k-l k-1 fc-1 

p— 1 a— p—1 a— 

Also, to go from Branson's formulas, we applied the functional equation 

Ci/2(s) = (2^ - l)Cfl(s), (7) 

since it turns out convenient to express everything in terms of the Riemann zeta 
function C_r. 
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We shall repeatedly rely on the controllability of (^r along the positive real 
axis. In fact we shall apply only the following two quite elementary facts. 

< CrH <Cr, ne N+, (8) 
0<CkH<C'j„ neNo, (9) 

for suitably chosen constants. 

If now n is odd, we have as always C{d^,S"){(^) — 0, since the kernel is trivial. 
Thus we let n be even and apply ([5]). Recalling the functional equation for the 
Riemann zeta function, in the form most convenient here, 

Cr{s) = 2^^^-i sin - s)Cr{1 - s), (10) 

we may now rewrite in a form where the asymptotics are more readily seen. 

C(i5^5-)(0) - E rfa,fe(-l)"+'(2^)-'"-'(2« + l)!Cfl(2a + 2). 

^ a=0 

Note that for fixed k this is a sum in which all the terms have the same sign. 
Namely 

sign = (-l)'=-i-", 

while (li is positive, since C-r(O) (= is excluded. One thing we get from this 
is 

sign [C(D2,s-)(0)] =(-1)'. 

Very importantly the constant sign (for fixed k) allows us to estimate term by 
term and we write 

|C(zj2,s-)(0)| <C^A(fc), (11) 
where the positive numbers A{k) are defined as 

m ■■= (^^^(-1)'+' E da,.(-l)"(27r)-2"-2(2a + 1)!. (12) 

The point is now that we can control the indirectly defined polynomial coeffi- 
cients da^k through a simple recursion. From ((6]) we get 

da,k+l — da-l,k — k^da.k- (13) 

Note that it is of course implicitly understood that da^k = 0, if fc is not in 
the range {0, 1, . . . , fc — 1}. By a change of index for the first term and using 
2Q!(2a + 1) < 2fc(2fc + 1), this shows the following estimate. 

fc 1 



2fc + 1 27r2 



A{k) (14) 



< lA(k). 
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Thus for instance with 5^ = | < 1, we get 

C(D^5")(0) = O(<5")> 

proving in particular the claim that 

lim C(r>^S")(0) = 0, 



(15) 



meaning that the phase converges to zero. Furthermore note that with = 
0.55 . . . this convergence is indeed rapid, as also witnessed by the numerics of 
Bar-Schopka. 

To deal with the derivative we differentiate (O at s = and use C(~2a) — 
0, a e N+. 



2^+2 



fc-i 



C('d2,s2'»)(0) = , ^rfa,fcCfl(-2a- 1), 



lk+1 



(0) = 



(2k - 2) 



eo.fe 



fe-l N 

log2 + ^e„,fc(2-2"-l)CA(-2a)l. 



Again we use the functional equation (|10p which by differentiation gives 

CA(-2a - 1) = 2(-l)"(27r)-2"-2(2a + l)!CA(2a + 2), a e Nq, 
+ [log(27r)+7-i72a+i]Ci?.(-2a-l), 

CA(-2a) - ^(-l)"(2^)-2"-i(2a)!Ci^(2a + 1), a G N+. 
Here 7 is Euler's constant, and we have applied 

n-l 



r{n) 



Hn-i-j, i?„_i = Vi, neN. 
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(16) 
(17) 

(18) 
(19) 

(20) 



Dealing first with the even dimensional case n = 2fc, we note that the second 
term in (|18p gives a contribution that converges to zero. This follows since in 
absolute value, after summing over a, it is altogether dominated by 



(log(2^) + 7 + H2k+i)CiD-^s-^){0) = OiS^"") 



where again 5^ = | < 1 for example is admissible, since {^(1 + l/7r^)}'°i/2*;+i 
is bounded. 

To control now the first term in ([IS]), note with A{k) from (jl2p we have similar 
estimates here 



(2fc- 1)! ^ 

^ ' a=0 



< C'j,A{k), 



thus we have the desired convergence 

lim C(D2,s2'=)(0) = 0, 
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and in fact C^di 526) (0) — 0((5^'^) with, for example, (5^ = | again. 

In odd dimensions the situation is changed only slightly. The polynomial 
coefficients are now ea^k from ([6]) and as before we have 

sign[e„,fe] = {-If-^'". (21) 

Writing now ([T7]) as 

C(D2,s2'=-i)(0) = - (2fc_2)! 51 ea,fc2'(-2a), 

with 

Z(-2.) (22) 
^ ^ |log2 a = 0. ^ ^ 

Note that once again the signs match up 

sign[Z(-2a)] = (-I)''-!-", 
and for a constant Cr :— max (log 2, ^) the estimates 

\Z{-2a)\ < C(27r)-2"(2a)! 
hold, and we define a new sequence of positive numbers B{k) by 

|C('z..,S--)(0)| < m ^|^(_l)'=-i|^'e„,,(-l)"(27r)-2"(2a)!. 

Again we find recursion relations for the relevant polynomials 

eQ,fe+i = -(A: - i)^eQ,fe + eQ_i^fc, (23) 
which gives recursive estimates on the B{k) 

i?(fc + l)<{i + ^}i?(fc)<^i3(fc), 

and this concludes the proof of the theorem. 



3. Laplace operator sphere determinants 
Following [Dowlj we investigate Laplacians of the form 

L^ = A+'^^R-al, (24) 

where R is the scalar curvature (of S*") and the q;„ e M are constants, but may 
depend on the dimension n. We restrict here to the case < q;„ < ^^^-j^- Recalling 
now that 

R{S:,J = n{n-l), 

we see that a„ = ^^^-^ gives the ordinary Laplacian, while a„ = ^ corresponds 
to the Yamabe operator in each dimension. Note that there is no kernel of the 
operator L in dimension n if a„ < , and that for — 2i=J- the kernel is the 
constant functions on S"". 

The theorem we will prove in this section is the following. 
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Theorem 2. On the standard (radius of unity) spheres we have the limits 

lim det(A^can+') = 0. (25) 

k — >-oo 



lim det(y,^^l,J = l (26) 

n— >oo 

In fact, we can display the asymptotics as (for n — 2k 1) 

Cl,,„,o).io,„ + o(i2i^), 



' log log n 

and (for any n ) 



Remark 2. The proof of the claimed limit for the conformal Laplacian is easily 
carried out analogously to the Dirac case. However we give a different type of 
proof, which works for the ordinary Laplacian as well, and is more suitable for 
exposing the significance of the coupling constant and kernel of the operator. 

From [Dowlj and [Dow2] we get the corresponding zeta functions and their 
first derivatives at zero 



Cl(0) =C(0, a_ - Q!„) + C,'i(0- a+ - "") + C«(0, a- + a„) + C(0, «+ + "«) 

(27) 
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In the notation here, the numbers a± := (n± l)/2 originate in the contributions 
from Dirichlet (a+) and Neumann (a_) boundary conditions on hemispheres. For 
the ordinary Laplacian, the four functions a± ± a„ are, by order of appearance 
in the above equation, (0, 1, n — 1, rt). For the Yamabe operator it reads (§ — 

1 H li 21 _i_ 1 A 
-•^1 2 ' 2 ' 2 

The notation "[+ ln(n— 1)]" means that this term is only present if Lq, has non- 
trivial kernel. The numbers N2j{n) are defined using the following polynomials, 
namely let 

n-2 n-2 
p—l r—O 

Then 

N2j{n) := j^^N2,-2in). 

Finally the zeta functions appearing here are (special cases of) Barnes zeta 
functions, where for a € M 

C„s,a-E 1 r— Res>n, 28 

■■^ V n—1 I (a -\- mr 

Note that here, if a = 0, it is implicit that the term m = is omitted. 
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The meromorphic continuation is carried out using the contour integrals 

C„(g,a) = ^-^(^^^) / ,--idz, (29) 

' ' 2tt Jh (1 -e^)" ^ ' 

where is a left Hankel contour, and is defined using the negative real 
axis branch cut of the logarithm enclosed by H. The contours will be chosen 
conveniently for getting estimates of the Barnes functions. Deformations of the 
curves must of course respect the branch cut and the possible poles at 27ri/c, k e 
Z. 

We note a few general features of the Barnes zeta functions. Each contour will 
be taken as the positively oriented boundary of some box (— oo,r„] x [— 
where r„ is a real- valued function of the dimension n. We denote such a contour 
by 7„ and decompose in the obvious way into one vertical and two horizontal 
components /i+ U U /i~ . Now, since 

|l-e^|" = (l + e2^-2e^cos(2/))^ (30) 
is the norm of the denominator at z = x + z?/, we always have 

|l-e"|" = (l + e^^)^ (31) 
on the horizontal parts of the contours. From (|29p we see, given that a < n, 



(1 - e^)" 



< {l + e^'^'') 2^0, for n~^oo, (32) 



for fixed z and fc, i.e. pointwise convergence to zero of the integrands. Lebesgue 
dominated convergence theorem can therefore in most case be used to prove 
that the contribution from a horizontal piece is 0{n~°°) as n ^ oo, meaning by 
definition 0{n~^) as n ^ oo, for any k. 

When a > 0, the derivatives at s = are, using ([20|l 

C„(0,a„) = — / — dz. (33) 



oo, 



27r y^,^ (1 ~ e^) 
Proposition 1. //a„ — 0, 1, ^^^^ or ^, then 

C^(0,a„)^0, as n 
in fact it is 0(n~°°) as n oo. 

Proof. For dealing with these cases we use the fixed contour consisting of the 
positively oriented boundary of the box (— oo,log4] x [— ^,^]. From the gen- 
eral remarks following ([5^ . the horizontal parts of the contours all contribute 
0{n~°°) for a„ > 0, by Lebesgue dominated convergence. 

In the case a„ = however, note that the continuation of the integral doesn't 
work directly around s = 0. From (j28|) we deduce, by straightforward resumma- 
tion, a family of functional equations 

Cn(s, 0) = Us, 1) + (n - l)C„(s + 1, 1), (34) 
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thus allowing again for meromorphic continuation to C by integral representa- 
tions. Taking into account the front factor of 7^(1 — s) in ([29|) and using 

r(i - s) = - 7 + . . . 

I — s 

near s — I, shows 

Cis, 0) = c (0, 1) + ^ jrhf ^ + ^ (35) 

On the vertical line segments u„ we use 1 1 — | > 3 and estimate 



(n-1) 



(1 - e^)" 



P{z) dz 



< Cp(n- 1) 3-" ^ 0, as n^oo, 



where P is any polynomial in - and logz, as in the expressions appearing in 
(1331) and (El). 



In the remaining cases a„ — -^^^ , we focus for the sake of argument on the 



2 ' 2 
latter, i.e. 



log z + 7 



jj (e-/2 „ z 

Then using | e^^^ — e~^/^|^ =e^-|-e^^ — 2 cos(j/) we see that indeed, by Lebesgue 
dominated convergence, there is rapid convergence to zero, namely 

C(V,S")(0) = O(2-"), 

concluding the proof. 

The next proposition deals with the most complicated of the terms, with param- 
eters a„ = n — 1 and n. Here there is no way to deform the Hankel contour in 
the complex plane, respecting the poles and branch cut of logz, so as to obtain 
pointwise convergence to zero of the integrand as n — s- everywhere along the 
curve. This follows from ()30p . since it would imply a; < at some point of the 
curve, which would thus intersect the logarithmic branch cut. 

The following lemma essentially gives a suitable hybrid of Laplace's method 
and the method of stationary phase, applicable to the cases needed here. 

Lemma 1 (Laplace's method/stationary phase hybrid). Assume ip : [a, b) - 
K>o for < a < b < 00 satisfies 

(i) if is continuous and bounded, 

(ii) 2; I— + • tf{x) is increasing on [a, 6). Then, writing sjj := tan~"'^(e~''), 



b 

{ie~^ — 1) (p{x)dx 



<— ^ max ^(log(^))+0(n-°°), 
Titansfc ^^<^<^^+2L 



where the error term is uniform in b, but not generally in a. 

Remark 3. (i) Note that in applying this lemma, we will let b depend on n, 
while a will be fixed. 

(ii) If is C^, the second condition in the lemma is equivalent to ip' > —ip. 
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Proposition 2. //a„ n — 1 or n, then 

C(0, a„) 0, for n ~> oo. 

Proof. For proving this, we assume n > 4 and shift for each n the contour to 
the positively oriented boundary of the box (— cx),logn] x [— which stiU 
encloses only the singularity at z = 0. 

The contributions from with, say, Rez < log 4, are again 0{n~°°) by 
while on the moving right hand edges w„, we have the estimates 



e 



ir"-(i + ;^-f-(^)r^<(i-^)""> 

yielding, for some constant C > 

,log logn 



1 log z + 7 

az 



< C- 



\ogn 



Thus the contribution from this part tends to zero as n — > c», and similarly for 
cin = n — I, again with or without the logz present, as needed for the terms in 

For the contributions from the right half-plane part of the horizontals, denoted 
by in the cases a„ = n—l,n, any finite piece gives an 0(n"°°) contribution 
as n ^ oo. 

To deal with the case a„ = n, we calculate explicitly, 

1 /-'"g" _ ,_J^- W2) ( log ^x^ + i^y + i tan-i (i:) + 7) 
Im / (ze "=-1) ^ 7— -r -dx. 

To take the oscillation into account, we repeatedly apply Lemma [T] to the ex- 
pression, with Sfc = i, considering the product terms separately. For each term 
the relevant positive function is decreasing for x large (i.e. for n large, since 
finite pieces contribute 0{n^°°)), so the maximum is evaluation at Sb + ^- For 
example one term is analyzed as follows 



xiogy^^+TfF 



'^i(^) ■= ^2 ^. (1)2 — ' a;e[i,iogn), 



log n 



so that by Lemma [Jl the righthand side converges to zero, as the remaining 
terms can similarly be shown to do. 

In the case a„ = n— 1, Lemma[l]is not needed, since after a partial integration 
it is easily seen that 

(e--l)""e-- ^°^^ + ^ rfz = 0(n-^) as n 00, 

HL + [0,logn] Z 

which ends the proof of Proposition [2l 
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Remark 4- Note that the use of Lemma[T]in the proof was essential, since without 
the oscillating factor, amounting approximately to sin(ne~^), we would have 

i(log(logn))^> T'" (l + e-2-)-"*'i^dx>i^{log(logn)-ilog2}, 
and thus convergence to +00, at a rate between log log n and log^ logn. 



We need to deal with the last term in ([27|) . which is controlled as follows. 
Proposition 3. // n is odd, then 

Setting a„ = ^ , then for n of any parity 

LfJ 3-1 



lini y ^iV2, (n) V 



1 



foemi? 0(2 ") as n ^ 00. 



Remark 5. Interestingly enough, this term is exactly the "correction term" from 
a product formula of certain determinants discussed in [Dowlj , and this makes 
it particularly interesting to see that this tends to zero. The convergence is 
exponential and this means that the anomaly, in this situation quickly evaporates 
as n — > 00. 

Proof. If n = 2A; + 1 is odd we exploit symmetry and rewrite the left-hand-side 
in the defining equation for Nr{n) 

k-l 2k-l 



p—1 r—0 

showing the vanishing of all terms with even r. 
If n = 2A: is even we again use symmetry to write 

p—1 r—0 

This shows that sign[iV2,._2(2A:)] = (-1)^+'- and gives as usual a recurrence 
relation 

N2r-2i2k + 2) = N2r~ii2k) ~{k- \fN2r-2{'2k). 

We estimate 



j=i i=0 



(2^-1)! 
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Letting q:„ = i we find the recursive estimates 

J=i ^ (36) 

for sufficiently large fc, proving exponential convergence to zero and ending the 
proof. 

4. Determinants on rescaled spheres 

It may be argued that taking the spheres with the radius unity standard metrics 
is not as natural for the problem of determinants, and that rescaling to unit 
volume is more appropriate. For the conformally covariant Yamabe and Dirac 
operators, this is of course produces no change, and focus remains on the ordinary 
Laplacian. Here the rescaling indeed cancels the leading term logn in the zeta 
derivative at zero. 

If we rescale the metric g by the constant A > to ^ := X^g the determinant 
of the ordinary Laplacian changes as 

(0) = C' (0) + 2 log A • C(0) - C (0) - 2 log A, (37) 

where the last equality is due to ^(0) = — 1. To rescale the spheres we use 



A(n) = (vol(5"))--. I j.(B±i) ; ' 

which with Stirling's formula 

V2^Ti"+^e""+T2^ < n! < \/2^n"+5e-"+ 

gives 
logA(r7,) 

lrn+1 -I, _ n, n— 1 n— 1 1 -i „, 



n 3) 



i logn- 1(1 + log(27r)) + 0(n-i). 



Inserting in (|37|) gives 

Ccscaicd(O) = C'(0) - logn + 1 + log(27r) + 0{n-'), 

thus cancelling the leading order term log n from the case of radius unity spheres. 
By the above along with Theorem [21 we find the new limit and asymptotics. 

Corollary 1. On the standard spheres rescaled to unit volume, we have the limit 

j^li^i^dct(Z\,5^^<!;+\od) ^ 
In fact the asymptotics can be displayed as (for n = 2k + 1) 

' log log n ^ 



Ck.^(0)^l + Iog(2.) + O(ifi^). 

V logn / 
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5. Proof of a stationary phase lemma 

Proof (Proof of Lemma[I\) . Let Sq :— ta.n~^ {e~°'), si, :— ta.n^^{e~^) and compute 



sm 



(ntan-i(e-^))(l + e-2^) ^ ip{x)dx 



sm 



(nx) (l + tan^ x 



-^^(log(dr^)) 



tanx 



dx 



E 



sm 



{nx) (l + tan^ x) 2 



-^^(log(tsb)) 



(38) 



sin (nx^ (l + tan^ 



tan a: 



tanx 



dx, 



where n, 



Sb + - 



b • = 



l + tan^x) 2 BVtan:.;; 



Thus the last term above can be estimated by 
dx < C^^a ( 1 + tan^ (sa ) 



tanx 



rt-2 
2 



when n > J- , and is 0(n~°°) as n — > cx), for any /c. Note as claimed, that the 
error term is only uniform in the parameter b. The first term in (|38p is written 
as a sum over the half-periods of the sine function. Fixing n, the alternating 
behavior gives the estimate 



E/ ffe"n, sin(nx)(l + tan2x) 



-^^(log(t5b)) 



fe=i 

■Sb + 



tanx 



dx 



< 



< 



— - — max 

tan Sb Sb<x<Sb+ 



in (nx) I (l + tan^ x) 

¥'(iog(dnj)) 



^^(i°g(dnj)) 



tanx 



dx 



Similar computations apply for the contribution from the real part of the integral. 
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